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ABSTRACT 

High  speed  Mercury  Cadmium  Telluride  photovoltaic  detectors  sensitive  to  infrared 
emission  are  investigated  as  a  means  of  measuring  real-time  surface  temperature  on  a 
micro-second  timeframe  pulsed  ablative  discharge  with  Teflon™  as  the  ablated  material. 
Analysis  is  used  to  derive  a  governing  equation  for  detector  output  voltage  for  materials 
with  wavelength  dependent  emissivity.  The  detector  output  voltage  is  experimentally 
calibrated  against  thermocouples  embedded  in  heated  Teflon.  Experimental  calibration  is 
performed  with  Teflon  that  has  been  exposed  to  ~200  pulsed  discharges  and  non  plasma- 
exposed  Teflon.  These  are  compared  to  theoretical  predictions  to  analyze  emissivity 
differences  between  plasma-exposed  and  non  plasma-exposed  Teflon.  The  diagnostic 
capability  is  evaluated  with  measurements  of  surface  temperature  from  the  Teflon 
propellant  of  electric  micro-Pulsed  Plasma  Thrusters.  During  the  capacitive  discharge, 
there  is  insufficient  information  to  claim  that  the  surface  temperature  is  accurately 
measured.  However,  immediately  following  the  discharge,  the  post-pulse  cooling  curve 
is  measured.  The  statistical  spread  of  post-pulse  surface  temperature  from  shot-to-shot, 
most  likely  due  to  arc  constriction  and  localization,  is  investigated  to  determine  an 
operational  envelope  for  post-pulse  temperature  and  mass  ablation.  This  information  is 
useful  for  determining  post-pulse  ablation  contributions  to  mass  loss  as  well  as  evaluation 
of  theoretical  discharge  models  currently  under  development. 
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1.  INTRODUCTION 


Propellant  surface  temperature  measurements  are  critical  to  further  development 
of  the  micro-Pulsed  Plasma  Thruster  (micro-PPT).1  A  micro-PPT  is  a  spacecraft 
propulsion  device  that  uses  a  pulsed  (~20  us)  electrical  surface  discharge  across  the  face 
of  a  solid  Teflon™  propellant.  The  discharge  ablates  a  small  amount  of  material  (~4  pg), 
ionizes  it  to  plasma,  and  electromagnetically  accelerates  it  to  create  thrust.  Following  the 
pulsed  discharge  another  mass  (~1.2  ug)  of  neutral  vapor  and  particulates  is  emitted  from 
the  device.  The  surface  temperature  measurements  are  needed  to  better  understand  the 
conversion  from  solid  propellant  to  plasma  accelerant  and  characterize  the  effects  of  late¬ 
time  neutral  vapor  generation  to  thruster  performance.  In  addition,  post-pulse 
temperature  measurements  are  needed  to  develop  and  verify  models  of  interaction 
between  the  spacecraft  and  the  exhaust. 

Analysis  of  the  post-pulse  mass  ejected  from  PPTs  has  typically  been  performed 
with  interferometry  sensitive  to  neutral  vapor  after  the  pulse.  With  the  development  of 
the  micro-PPT,  the  low  energies  typical  of  the  discharge  have  limited  the  neutral  density 
to  levels  below  the  sensitivity  threshold  of  laser  interferometers.  A  Herriott  Cell  was 
used  to  augment  the  resolution  of  these  interferometers,  but  post-pulse  neutral  densities 
were  still  unmeasurable.  An  alternative  method  of  investigation  is  proposed  to  attempt  a 
direct  measurement  of  the  propellant  surface,  since  the  equilibrium  vapor  pressure  curve 
for  Teflon  is  well  known.1,4 
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Direct  measurements  of  the  propellant  surface  temperature  can  be  achieved  using 
photovoltaic  cells.  Research  in  the  fracture  dynamics  community  relies  heavily  on 
infrared  photovoltaic  detectors  to  characterize  temperature  rise  near  the  tips  of  fast- 
propagating  cracks  in  solid  materials."’  ’  The  time-frame  under  observation  for  these 
experiments  is  similar  to  timeframes  during  which  a  micro-PPT  discharge  is  observed.  It 
is  a  claim  of  this  paper  that  real-time  propellant  surface  temperature  measurements  can  be 
accomplished  using  the  similar  techniques.  The  relationship  between  surface  temperature 
and  Teflon  vapor  pressure  is  well  known  at  high  temperatures.  Through  this  relationship, 
surface  temperature  (Ts)  governs  ablation  rate  in  these  thrusters  and  becomes  critical  in 
calculation  of  total  mass  loss,  available  charge  carriers,  and  late-time  ablation.  These 
extend  the  influence  of  Ts  to  plasma  resistivity,  ionization  fraction,  current  density,  and 
the  feedback  process  of  ablation  from  arc  radiation  and  particle  convection.  Surface 
temperature  is  the  defining  parameter  for  thruster  operation  and  analysis,  which  until  now 
has  been  un-measurable  in  real-time  for  these  thrusters. 

This  paper  is  directed  towards  validating  the  use  of  the  photovoltaic  detectors  in 
high-speed,  pulsed  plasma-wall  interaction  applications.  An  inherent  problem  for 
performing  this  measurement  is  the  wavelength  dependent  emissivity  of  polymers  used 
for  pulsed  discharge  applications.  Teflon  has  a  wavelength  dependent  emissivity  and  is 
used  for  a  number  of  applications  including  electric  rocket  propulsion  ’  ’  and  plasma 
opening  switches.10’11  Teflon  is  also  an  ideal  candidate  material  for  this  investigation. 
Vaporization  from  the  surface  of  Teflon  generally  leaves  a  smooth  surface  with  minimal 
changes  in  emissivity. 
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Section  2  describes  the  apparatus  and  experimental  approach  for  calibration  of 
photodiode  voltages  from  heated  Teflon  along  with  diagnostic  application  to  the  thruster 
face.  Emissive  differences  between  plasma-exposed  and  non-plasma-exposed  Teflon  are 
addressed  experimentally.  In  Section  3,  a  theoretical  estimate  of  the  expected  detector 
output  for  heated  Teflon  is  derived.  Data  from  infrared  spectroscopy  is  used  to  estimate 
wavelength  dependent  emissivity.  Theoretical  prediction  helps  assess  background 
emissivity  differences  between  plasma-exposed  and  non-plasma-exposed  Teflon  samples 
for  calibration.  Additionally,  a  short  review  of  the  Keidar-Boyd  Model12"14  predictions 
for  surface  temperature  is  presented.  In  Section  4,  the  diagnostic  capability  is 
demonstrated  with  measurements  of  surface  temperature,  comparison  with  Keidar-Boyd 
Model  predictions,  and  discussion  of  the  limitations  of  the  diagnostic. 

2.  APPARATUS 


2,1  Plasma  Source 

The  plasma  source  is  a  micro-Pulsed  Plasma  Thruster1  with  an  outer  electrode 
diameter  of  6.35  mm,  an  inner  electrode  diameter  of  1.64  mm,  an  outer  Teflon  diameter 
of  5.46  mm,  and  RLC  discharge  energy  of  4.35  J.  The  source  is  operated  in  single  pulse 
mode  to  minimize  steady-state  temperature  gain  throughout  the  thruster  body  and 
propellant. 
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2.2  Detector  and  Optics 


The  HgCdTe  photovoltaic  detector  is  chosen  for  its  broad  range  of  responsivity 
over  2-12  pm  wavelengths  and  high  speed  response.  A  linear  array  of  4  detectors  made 
by  Fermionics15  has  DC  to  10  MHz  bandwidth  for  an  active  area  of  80  pm  x  80  pm 
square  and  a  pitch  of  20  pm  between  detectors.  Peak  responsivity  (Rp)  is  8.2  A/W  with  a 
detector  rise  time  ~40  ns.  Figure  1  shows  the  4  detector  array  superimposed  over  the 
thruster  face  to  scale  for  assessment  of  the  spatial  resolution. 

The  optics  for  infrared  detection  use  a  15.2  cm  diameter  gold  coated  concave 
spherical  mirror  with  a  30.5  cm  focal  length  and  a  50  mm  diameter  gold  coated  flat 
mirror.  A  distance  of  61  cm  from  the  spherical  mirror  to  both  the  target  and  the  detector 
is  maintained.  The  50  mm  flat  mirror  is  placed  36  cm  from  the  spherical  mirror  to 
complete  the  Newtonian  telescope.  Figure  2  shows  a  schematic  of  this  layout. 

Targeting  for  both  the  calibration  apparatus  and  thruster  face  during  experimental 
testing  is  performed  using  crossed  tungsten  wires  of  80  pm  diameter.  Current  up  to  1 .2  A 
is  run  through  each  wire  to  induce  thermal  emission.  The  signal  is  maximized  by  moving 
the  detector  array  forward  and  backward  to  locate  the  signal  maximum  corresponding  to 
the  focal  point  of  the  telescope.  The  Newtonian  telescope  is  then  used  to  target  first 
horizontally  and  then  vertically  with  the  two  crossed  wires. 
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2.3  Calibration  Assembly 


To  calibrate  detector  voltage  output  with  actual  heated  Teflon  temperature,  a 
3.175  mm  thick  sheet  of  Teflon  is  sandwiched  between  two  copper  plates  with  optical 
access  holes  5.38  mm  in  diameter  in  both  copper  plates.  Three  type  K  thermocouples  are 
drilled  into  the  sides  of  the  Teflon  sheet  and  embedded  close  to  the  point  where  the  IR 
detector  is  imaged.  The  Newtonian  telescope  images  the  detectors  at  1 : 1  magnification 
onto  the  Teflon  surface.  Figure  3  shows  a  photo  of  the  front  of  the  calibration  apparatus 
with  thermocouples  embedded  (left)  and  a  photo  of  the  apparatus  with  the  back  removed 
highlighting  the  thermocouple  shadows  embedded  in  the  Teflon  sample  (right).  This 
illustrates  the  relative  location  of  the  thermocouple  measurement  and  the  detector  focal 
point.  The  baseplate  of  the  apparatus  is  seated  on  a  Micropyretics  Heaters  International 
(MHI)  hotplate.  Two  copper  plates  are  hinged  to  form  a  triangle  of  reflective  metal  that 
sits  over  the  apparatus  and  reflects  heat  back  at  the  calibration  apparatus.  Optical  access 
for  the  Newtonian  telescope  is  provided  through  a  25.4  mm  conflat  viewport  made  of 
Zinc  Selenide  (ZnSe).  Heat  is  transferred  through  the  baseplate  to  the  front  copper  plate 
which  then  heats  the  front  side  of  the  Teflon.  Data  from  the  thermocouples  and  detectors 
are  recorded  by  a  data  logger  at  1  or  2  second  intervals.  A  typical  calibration  takes  about 
10  hours.  A  slow  heating  scheme  is  used  to  increase  accuracy  by  reducing  transient 
thennal  effects. 

During  pulsed  discharge,  problems  with  Teflon  surface  charring  "  can  occur  that 
affect  material  emissivity.  There  also  exists  the  possibility  that  a  clean,  plasma-exposed 
surface  will  exhibit  different  emissive  characteristics  than  a  virgin  (non  plasma-exposed) 


6 


Teflon  surface.  A  demonstration  of  the  effects  of  the  plasma  on  a  plasma-exposed 
surface  is  necessary  to  detennine  if  the  exposure  affects  the  emissivity  of  Teflon 
propellant,  and  therefore  the  calibration  curve.  A  thruster  breadboard  with  copper 
electrodes  and  a  Teflon  sheet  underneath  is  used  to  expose  Teflon  suitable  for  calibration 
tests  to  plasma  conditions  experienced  by  the  micro-PPT  surface  to  determine  if  there  are 
significant  changes  to  IR  emissivity.  Samples  have  been  exposed  to  4.35  J  discharges  for 
pulse  numbers  ranging  from  100  up  to  1500  pulse  exposures.  These  are  compared  with 
non  plasma-exposed  Teflon  samples  to  characterize  changes  in  IR  emissivity  that  may 
occur  with  plasma  exposure. 

2.4  Quantifying  Plasma  Effects 

To  experimentally  characterize  the  emission  from  the  plasma  and  neutral  vapor 
above  the  fuel  face,  IR  measurements  are  taken  with  the  detector  optical  cone  targeted  to 
a  point  in  front  of  the  central  electrode  and  traversing  across  the  exit  plane  of  the  thruster. 
The  view  cone  in  this  orientation  does  not  intersect  the  surface  of  the  micro-PPT.  Figure 
4a  shows  the  cone  apex  targeted  schematically  to  illustrate  this  configuration.  Figure  4b 
shows  the  apex  of  the  optical  cone  intersecting  the  propellant  face,  the  configuration  used 
when  measuring  the  surface  temperature  of  the  propellant.  By  measuring  emission  from 
the  plasma  and  neutral  vapors  without  a  heated  background,  the  emissive  contribution  of 
the  plasma  and  neutral  vapor  is  characterized.  If  the  plume  contribution  ends  when  the 
current  pulse  ends,  then  the  surface  can  be  measured  after  that  time  without  emissive 
interference  from  the  exhaust.  To  characterize  this,  measurements  across  the  fuel  face 
are  made  with  the  focal  point  of  the  optics  located  at  1  mm,  3  mm,  5  mm,  and  7  mm  in 
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from  of  the  fuel  face.  When  these  measurements  are  compared  to  fuel  face 
measurements,  an  assessment  concerning  optimal  viewing  time  during  the  pulse  is  made. 
Repeated  tests  show  that  after  the  current  reversal  stops,  there  is  no  measurable  detector 
signal.  These  results  are  discussed  below,  but  here  it  is  sufficient  to  state  that  after  the 
current  pulse  ends,  plasma  and  vapor  emission  contribution  also  ends  and  a  clear  surface 
measurement  can  be  obtained. 


3.  THEORY 


3.1  Teflon  Emission 

The  distribution  of  radiation  emitted  from  a  body  depends  on  temperature  and  is 
described  by  the  Plank  function16 


PE  = 


Q 


A5 


(  cy 

yAT 


-1 


(1) 


where  Cl  is  3.742xl08  W-pm4/m2,  C2  is  1.439xl04  pm-K,  A  is  the  radiation  wavelength 
(pm),  and  T  is  the  material  temperature  (K).  This  equation  yields  spectral  emissive 
power  for  a  black  body.  For  Teflon  emission  we  derive  an  expected  detector  output  from 
calibration  with  a  Teflon  sample  using  data  from  transmissive  spectroscopy  as  a  starting 
estimate.17  For  any  material,18 

p  +  a  +  x  =  1  (2) 

where  p  is  the  reflectivity,  a  is  the  absorptivity,  and  x  is  the  transmissivity.  Areas  of  high 
absorption  correspond  to  areas  of  high  emission  since  energy  absorption  equals  energy 
emission  from  the  stretching  and  bending  modes  of  the  C-F  bonds  in  the  Teflon 
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molecules.  The  reflectivity  of  Teflon  in  the  wavelengths  of  interest  is  unknown,  but  from 
Eqn.  2  it  is  assumed  that  material  reflectivity  decreases  the  amount  of  total  emission. 
Therefore,  if  the  reflectivity  is  omitted  in  a  prediction  of  detector  response,  an 
overprediction  of  the  total  emission  is  expected.  This  expected  overprediction  in  noted  in 
the  final  estimates.  Eliminating  reflectivity  and  setting  absorptivity  equal  to  emissivity 
(s)  converts  Eqn.  2  to 


8  =  l-i  (3) 

allowing  an  estimate  of  the  wavelength  dependent  emissivity  for  Teflon.  An  emissive 
peak  of  unity  at  8.4  pm  is  the  chief  emissive  source4,17  along  with  an  overtone  at  4.4  mm 
of  ~0.5  due  to  superposition  of  the  stretching  and  bending  modes  present  in  the  C-F 
bonds.  Figure  5a  shows  an  estimate  of  Teflon  emissivity.  Using  this  estimate,  an 
expected  spectral  emissive  power  for  the  2-12  pm  range  is  calculated. 

A  first  order  detennination  of  the  expected  total  emission  and  response  is  obtained 
by  combining  the  Plank  distribution  for  Teflon  with  the  wavelength  dependent  detector 
responsivity  and  integrating  over  the  range  of  wavelengths  (2-12  pm).  This  gives  an 
expected  detector  current  output  calculated  by 


1  r.  =  J/DRpS>Prca  (4) 

where  D*  is  the  normalized  wavelength  dependent  responsivity,  Rp  is  the  peak 
responsivity,  s>„  is  the  wavelength  dependent  emissivity,  and  Pe  is  the  spectral  emissive 
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power  calculated  by  Eqn.  1 .  Equation  1  yields  the  total  power  emitted  from  the  specimen, 
however,  only  a  fraction  of  that  power  is  observed  by  the  detector.  Assuming  1:1 
imaging,  with  a  detector  view  angle  limited  by  the  flat  mirror  (50.8  mm  DIA)  the  half 
angle  of  the  viewable  cone  is  3.05°  on  average.  The  angular  area  of  a  hemisphere  is  2n, 
and  the  angular  area  of  a  small  circle  of  angular  radius  0  is  given  by19 

Clc  =  2;r(l  -  cos  0)  (5) 

The  ratio  of  the  circle  solid  angle  to  the  hemispherical  solid  angle  is  simply  (l-cos0), 
which  is  the  ratio  of  the  power  received  to  the  power  emitted. 

The  effect  of  the  optics  on  total  intensity  is  included.  The  Newtonian  telescope 
uses  two  gold  mirrors  with  a  reflectivity  of  0.99  in  the  relevant  wavelength  range. 
Additionally,  since  the  calibration  is  performed  at  vacuum,  the  transmissivity  of  the 
viewport  material  is  included.  The  optical  viewport  on  the  vacuum  chamber  is  ZnSe, 
chosen  to  maximize  the  total  transmission  in  the  LWIR  (t=0.7  from  2-12  mm).  The 
optical  response  is  defined  as  a  constant,  k0,  which  is  the  multiplied  value  of  all  mirror 
reflectivities  and  window  transmissivities.  Multiplication  of  the  ultimate  detector  output 
by  this  constant  corrects  for  losses  in  the  optical  system. 

Since  the  optics  provide  an  imaging  system,  the  detector  is  operated  with  a  1 : 1 
magnification  of  the  detector  element  onto  the  surface  of  the  specimen.  Only  emission 
coming  from  within  the  viewable  cone  (the  apex  of  which  is  the  detector  image)  can  be 
sensed.  The  detector  active  area  limits  the  total  number  of  photons  collected  from  the 
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optics,  therefore  the  active  area,  A,  is  included.  Each  detector  is  a  square  with  80  pm  to  a 
side  for  an  area  of  6.4x10'  m. 

Finally,  the  amplifier  has  a  transimpedance  factor  (Vout/Iin),  of  20000  V/A.  This  is 
the  gain  of  the  pre-amplifier  and  is  denoted  here  as  G.  The  combined  end  result  for 
estimating  detector  voltage  out  is 

V  =  koGA(l-cos0)jDtRp8,PedX  (6) 

x 

Using  Eqn.  6  and  noting  the  omission  of  Teflon  reflectivity  in  the  estimation  of  the 
emissivity,  a  calculation  for  the  output  voltage  from  this  detector  measuring  heated 
Teflon  can  be  made.  Since  an  overprediction  is  expected  for  the  Teflon  response,  an 
experimental  calibration  curve  is  obtained  to  demonstrate  the  deleterious  effects  of  the 
assumptions  up  to  this  point.  These  results  are  discussed  in  section  4. 

Comparison  of  the  predictions  from  Eqn.  6  with  experimental  calibration  requires 
the  estimate  of  the  emissivity  shown  in  Figure  5a.  To  detennine  if  calibration  results  are 
valid,  the  heating  and  cooling  curves  for  detector  voltage  out  vs.  surface  temperature 
must  match.  If  the  thennocouples  are  measuring  close  to  the  location  of  the  emission 
measurement,  and  there  are  no  changes  to  surface  condition  during  the  test,  both  curves 
will  demonstrate  the  same  emission.  Reasons  why  the  two  curves  might  not  match 
include  Teflon  expansion  out  of  the  viewing  hole  (which  causes  a  change  in  the  optical 
path),  targeting  errors,  material  depositing  on  the  surface  of  the  Teflon  causing  a  change 
in  emissivity,  or  possibly  a  permanent  change  in  emissivity  due  to  thermal  cycling  of  the 
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Teflon.  For  the  experimental  setup  used  here,  the  heating  and  cooling  curves  match  and 
are  used  to  determine  a  useful  estimate  of  emissivity.  With  all  other  terms  in  Eqn.  6 
known  as  well  as  the  major  emissive  structures,  the  emissivity  background  is  adjusted  to 
match  the  experimental  calibration  curve  discussed  in  Section  4. 1 .  This  allows  calculated 
extrapolation  of  temperatures  beyond  the  physical  calibration  limit  of  750  K. 

3.2  Vapor  Optical  Thickness 

The  experimental  setup  used  here  cannot  measure  surface  temperature  during 
current  reversal  due  to  emission  encountered  from  plasma  and  neutral  vapor  situated 
above  the  Teflon  face.  As  the  surface  is  heated  during  the  arc,  these  two  regions  exist  at 
high  temperatures  and  may  be  capable  of  emission  within  the  measurement  range  of  the 
detectors.  Determination  of  whether  surface  emission  may  be  viewed  is  dependent  upon 
the  optical  thickness  of  the  plasma  and  neutral  vapor.  Calculation  of  wavelength 
dependent  optical  thickness  is  beyond  the  scope  of  this  inquiry,  but  a  more  simple  optical 
thickness  calculation  is  possible. 

Presence  of  plasma  may  hinder  photon  collection  by  re-absorption  of  the  radiated 
photons.  The  magnitude  of  this  effect  is  estimated  by  calculating  the  optical  thickness  as 
the  non-dimensional  ratio  of  the  plasma  scale  length  to  the  characteristic  length  of  photon 
motion  prior  to  reabsorption. 
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Therefore  an  optical  depth  much  less  than  1  indicates  a  high  likelihood  for  photons 
generated  within  the  plasma  to  reach  the  edge  without  reabsorption  and  be  available  for 
collection  by  the  diagnostic  optics. 

The  plasma  optical  thickness  for  a  doppler  broadened  line  is  given  by  (CGS 

units) 

1 

( Mr2  V 

r  =  1.76x  10“13k  — —  N  /  (8) 

l  kT  J  P 

where  M  is  the  atomic  mass  (averaged  to  16.7  AMU  or  2.77x10'^  gm),  c  is  the  speed  of 
light,  k  is  the  Boltzmann  constant,  Te  is  the  electron  temperature,  Ne  is  the  number 
density,  and  lp  is  the  plasma  depth.  Inserting  typical  values  for  a  micro-PPT  (Ne=lxl016 

-5 

cm'  ,  Te  =  11500  K  (1  eV),  /  =  0.635  cm  (using  the  thruster  diameter),  6  pm  wavelength), 
gives  i  =  8xl04.  This  suggests  that  the  plasma  is  optically  thick.  Emission  from  the 
Teflon  surface  is  likely  to  be  reabsorbed  by  the  plasma  and  not  reach  the  detector.  The 
optical  thickness  of  the  plasma  may  limit  the  usefulness  of  the  diagnostic  to  after  the 
conclusion  of  the  plasma  discharge.  Note  that  this  analysis  for  optical  depth  presumes 
that  the  transitions  in  the  solid  Teflon  at  the  propellant  face  will  be  very  close  to  the  line 
transitions  in  the  plasma.  If  these  transitions  are  at  significantly  different  wavelengths, 
the  plasma  may  be  more  optically  thin  than  indicated  by  these  calculations. 

There  is  an  experimental  means  of  determining  whether  the  surface  emission  can 
be  independently  collected  during  and  after  the  pulse.  Both  plasma  and  neutral  vapor 
emission  can  be  accounted  for  experimentally  by  measuring  the  total  emission  in  the 
viewable  cone  in  an  area  where  they  are  present.  By  targeting  across  the  surface  of  the 
thruster  over  the  central  electrode,  the  plasma  and  vapor  contribution  can  be  assessed. 
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This  method  characterizes  the  gaseous  emission  throughout  the  pulse  and  is  described  in 
Section  2.4. 

3.4  Keidar-Boyd  Model 

The  Keidar-Boyd  model  is  an  end-to-end  (plasma  generation  to  plume  far  field)  model 
specifically  addressing  the  geometry  of  the  microPPT  with  regards  to  thruster 
performance  and  spacecraft  contamination  issues.  Three  models  are  employed  to 
investigate  the  complete  discharge  of  the  micro-PPT.  The  first  model  is  a  first-order 
ablation  model  that  characterizes  Teflon  loss  to  the  arc.  The  second  is  a  description  of 
the  Teflon  transition  through  neutral  vapor  to  plasma  in  several  sublayers.  Finally,  there 
is  a  Particle-in-Cell  (PIC)  description  of  the  plume  to  far  field.  These  models  use  a  fluid 
description  of  the  electrons,  a  direct  simulation  Monte  Carlo  (DSMC)  model  of  the 
neutrals,  and  a  PIC  description  of  the  ions.  The  ablation  and  plasma  generation  models 
provide  boundary  conditions  for  the  plasma  plume  model.12'14  This  allows  for  a 
consistent  description  of  the  plasma  flow  from  the  Teflon  surface  into  the  near  plume. 
Because  of  this  fact,  the  plume  properties  depend  very  heavily  on  the  upstream  conditions 
including  surface  temperature.  The  Keidar-Boyd  model  considers  arc  constriction  in 
calculation  of  surface  temperature  by  two  means  of  heating.  When  the  area  being  heated 
is  directly  under  the  constriction,  heating  is  performed  by  a  mixture  of  particle  convection 
and  radiation  to  the  surface.  In  cases  where  the  area  being  heated  is  far  from  the 
constriction  (i.e.  opposite  side  of  the  thruster  face)  the  dominant  heating  mechanism  is 
radiation  while  particle  convection  is  neglected.  Predictions  of  surface  temperature  from 
this  model  are  compared  to  experimental  results  below. 
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4.  RESULTS  AND  DISCUSSION 


4.1  Calibration  Results 

Figure  6  shows  the  calibration  data  for  both  non  plasma-exposed  and  plasma-exposed 
Teflon  along  with  the  theoretical  detector  output  using  the  emissivity  estimate  from  Eqn. 
6.  The  most  significant  result  is  that  the  small  degree  of  disagreement  between  virgin 
Teflon  and  plasma  exposed  Teflon  can  be  explained  completely  through  a  small  change 
in  background  emissivity.  Figure  5  a.  and  b.  show  the  difference  in  background 
emissivity  required  to  match  the  observed  difference  in  calibration.  For  the  non  plasma- 
exposed  case,  the  background  emissivity  is  set  at  0.04.  Keeping  the  major  emissive 
structures  constant  and  reducing  the  background  emissivity  to  0.01  results  in  good 
agreement  between  the  experimental  data  for  plasma  exposed  Teflon  and  the  modeled 
calibration.  Because  the  top  layers  of  Teflon  in  the  thruster  undergo  thermal  cycling 
from  repetitive  pulsing,  changes  in  the  crystalline  structure  result  in  a  slightly  different 
emissivity  than  non  plasma-exposed  Teflon.  The  small  degree  of  disagreement  between 
virgin  and  plasma  exposed  Teflon  calibrations  suggests  that  there  is  very  little  emissive 
contribution  that  cannot  be  traced  to  the  C-F  stretching  mode. 

Representative  data  points  from  experimental  calibrations  are  shown  along  with  a 
parabolic  fit  of  those  data  points.  These  are  also  compared  with  a  theoretical  prediction 
from  Eqn  6  (which  lies  directly  on  the  parabolic  fit).  The  fit  for  the  plasma-exposed 
calibration  curve  is 

V  =  2.87  *1CT5T2  (9) 
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Beyond  the  peak  measured  temperature  of  750  K,  the  voltage  is  calculated  using  Eqn.  6. 
For  these  experiments,  the  calibration  sample  is  conditioned  with  ~  100  pulses  before 
testing.  The  plasma-exposed  calibration  curve  is  used  to  translate  detector  voltages  into 
the  temperature  domain  for  measurements  from  thruster  propellant  faces. 

The  experimental  calibration  limit  of  750  K  is  due  to  solid  Teflon  deformation  at 
higher  temperatures.  This  has  the  effect  of  changing  the  optical  path  and  invalidating  the 
calibration.  Therefore  calibration  data  above  this  temperature  is  calculated  using  Eqn.  6. 
While  it  is  possible  that  the  calibration  could  change  at  these  temperatures,  there  is  no 
evidence  to  suggest  a  change  in  the  major  emissive  structures  at  8.4  and  4.4  pm  since 
these  are  the  result  of  vibrational  modes  present  in  the  polytetrafluoroethylene  monomers 
even  after  the  polymer  has  completely  unzipped. 

The  uncertainties  in  Figure  6  are  calculated  by  a  root-mean  square  standard  deviation 
of  5  calibration  tests  combined  with  the  ±0.3  mV  resolution  of  the  detector.  Also  shown 
are  temperature  uncertainties  using  the  same  method  based  on  the  thennocouple 
uncertainty.  For  the  K-type  thermocouples  used,  the  uncertainty  is  ±2.2  K  or  .0075 *T, 
whichever  is  larger.  Another  tenn  is  introduced  for  ±8  K  representing  the  temperature 
difference  between  the  embedded  thennocouple  location  and  the  center  of  the  calibration 
target  where  emission  measurements  are  made.  This  value  is  obtained  from  experimental 
data  and  ANSYS  predictions. 

4.2  Plume  Emissive  Contribution 

Measurements  of  the  plume  emission  alone  determine  the  contribution  of  plasma 
emission  during  the  discharge  to  the  total  signal  measured.  By  noting  where  in  time  the 
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plume  contribution  ends,  a  timeframe  for  possible  surface  viewing  is  determined. 
Measurements  taken  at  1,  3,  5,  and  7  mm  away  from  the  propellant  face  show  that  plume 
contribution  drops  to  zero  throughout  the  discharge  by  5  mm  from  the  fuel  face  indicating 
that  the  plume  expansion  at  that  point  has  limited  the  total  emission  that  the  detectors  can 
sense.  Figure  7  shows  an  uncalibrated  comparison  of  surface  measurement  with  a  plume 
measurement  at  1  mm  distance  from  the  propellant  surface.  Cessation  of  signal  from  the 
plume  occurs  at  or  before  the  end  of  the  current  pulse.  From  this  data  we  define 
acceptable  surface  temperature  data  as  that  which  occurs  after  the  current  pulse  is 
extinguished.  In  the  data  presented  below,  surface  temperature  data  begins  at  26  ps  after 
initiation  of  the  discharge. 

4.3  Surface  Temperature  Measurements 

Figure  8  shows  14  distinct  data  traces  taken  while  targeting  the  thruster  propellant 
face.  Initial  post-pulse  temperatures  (t=26  ps)  observed  range  between  600  and  850  K. 
Cooling  rates  appear  roughly  constant  from  pulse  to  pulse.  These  data  have  been  filtered 
digitally  using  a  low-pass  filter  cutting  off  frequencies  above  300  kHz.  There  is  a  small 
amount  of  electrical  noise  at  about  55  ps  not  removed  by  the  filter  present  in  each  pulse. 
This  defines  an  operating  envelope  of  surface  temperatures  that  can  be  explained  by 
differences  in  arc  spoking  characteristics  between  pulses.14 

The  spatial  resolution  of  the  detector  in  these  experiments  is  limited  to  the 
detector  size  as  imaged  on  the  propellant  face.  Since  this  is  1:1,  the  detector  only  sees  80 
x  80  pm  of  surface  area.  This  is  a  small  fraction  of  the  total  propellant  area  exposed  to 
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the  effects  of  the  arc.  Figure  9  shows  DICAM  images  of  a  micro-PPT  firing  at  a.)  2.3  J, 
b.)  5.55  J,  and  c.)  6.73  J.  These  firings  used  a  0.417  |uF  capacitor  with  the  central 
electrode  as  the  cathode.  The  discharge  voltages  are  then  a.)  3300  V,  b.)  5200  V,  and  c.) 
5700  V.  The  trend  observed  is  that  at  low  discharge  voltages,  there  is  a  greater  likelihood 
of  a  symmetrical  discharge  across  the  propellant  face.  At  the  higher  energies  the 
emission  suggests  that  anode  spots  form  on  the  outer  electrode  which  indicates  the 
possibility  of  current  constriction  across  portions  of  the  propellant  face.  If  the  detector 
measures  directly  underneath  an  arc  constriction,  the  post-pulse  surface  temperature  is 
expected  to  be  higher  than  is  the  detector  measures  a  portion  of  the  propellant  that  is 
distant  from  the  arc  constriction.  Using  the  voltage  and  temperature  uncertainties 
discussed  above,  individual  pulse  uncertainties  are  compared  for  two  of  the  extreme  cases 
in  post-pulse  temperature  measured.  Figure  10  shows  comparison  of  a.)  the  highest 
temperature  cooling  curve  from  the  data  set  (Pulse  20)  compared  with  the  surface 
temperature  predicted  for  an  area  covered  by  the  arc  constriction.  Part  b.)  shows  a 
comparison  of  the  lowest  temperature  cooling  curve  from  the  data  set  (Pulse  19) 
compared  with  the  surface  temperature  predicted  for  an  area  far  from  the  arc  constriction. 
The  Keidar  Boyd  model  predicts  that  a  temperature  measurement  made  in  a  location 
directly  under  the  constriction  will  have  the  highest  temperature  post-pulse  cooling  curve 
and  that  one  far  from  the  constriction  (i.e.  heated  only  by  radiation)  will  have  the  lowest 
temperature  cooling  curves  observed.  Any  spread  between  these  values  may  reflect 
changes  in  proximity  to  the  actual  constriction  location.  The  difference  in  these 
measurements  suggests  that  the  arc  was  moving  randomly  around  the  propellant  face 
during  testing.  This  is  not  conclusive  proof,  as  two-dimensional  imaging  of  the 
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propellant  face  in  the  long  wavelength  IR  would  be  required  to  assess  the  constriction 
location  positively  on  a  pulse-to-pulse  basis. 

The  natural  vapor  pressure  of  heated  Teflon  is  given  by  Wentink.4  The  mathematical 
description  of  this  curve  is1 

Peq=Pce^s  (10) 

where  peq  is  the  equilibrium  vapor  pressure  of  Teflon,  pc  and  Tc  are  characteristic 
pressure  and  temperature  values  respectively,  and  Ts  is  the  surface  temperature.  In  this 
case,  pc  =  1.84xl015N/m2  and  Tc  =  20,815  K.1  Combining  this  estimate  of  the  vapor 
pressure  with  the  measured  surface  temperature  yields  a  time-dependent  pressure 
calculation  relevant  to  neutral  vapor  generation  after  the  pulse.  Since  the  pressure  is 
logarithmically  dependent  on  the  surface  temperature,  small  changes  in  temperature  can 
result  in  large  differences  in  pressure.  Figure  1 1  shows  the  timeframe  from  26  to  40  ps 
expanded  with  post-pulse  pressure  shown  using  the  temperature  from  Pulse  20.  The 
pressure  uncertainty  is  based  only  on  the  uncertainty  in  the  temperature  measurement. 

Assuming  that  after  the  current  pulse  each  particle  liberated  is  a  CF2  monomer,  the 
average  particle  mass  is  50  AMU.  Multiplying  this  value  by  Avogadro’s  number  gives 
an  individual  molecule  mass  of  8.3x10"'  grams/particle.  Using  calculations  for  thennal 
evaporation  common  to  particle  vapor  deposition,  the  surface  temperature  and  vapor 
pressure  can  be  related  to  a  molecular  flux  through 

F=  ,  Peq  =  (11) 

yjh rmkTs 
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where  m  is  the  molecule  mass,  k  is  the  Boltzmann  constant,  peq  is  the  equilibrium  vapor 
pressure  (from  Eqn.  3),  and  Ts  is  the  surface  temperature.  This  relationship  is  derived 
from  the  Clausius-Clapeyron  equation.  The  molecular  flux,  F,  is  the  number  of 
molecules  liberated  from  the  surface  per  unit  area  per  unit  time.  Using  the  data  from 
Pulse  20  and  assuming  that  the  surface  temperature  measured  is  representative  of  the  total 
surface  temperature,  the  molecular  flux  can  be  multiplied  by  the  propellant  surface  area 
and  integrated  through  time  to  yield  a  total  molecular  evaporation  and  therefore  a  total 
post-pulse  mass  loss.  Figure  12  shows  the  integration  of  the  mass  loss  (and  molecular 
evaporation)  through  time  with  typical  uncertainties. 

The  total  mass  lost  through  post-pulse  evaporation  using  the  temperature  from 
Pulse  20  is  1.2±0.6  pg.  The  average  total  mass  lost  per  pulse  (5.2  pg/pulse)  is  measured 
by  pulsing  10000  times  and  measuring  thruster  mass  before  and  after  the  pulsing.  The 
post  pulse  mass  loss  as  a  percentage  of  the  total  mass  lost  is  23±1 1%.  This  assumes  a 
symmetrical  current  sheet  with  exposed  Teflon  propellant  area  of  2. 13x10"  m  .  Note  that 
the  assumption  of  the  measured  surface  temperature  across  the  propellant  face  does  not 
necessarily  apply  due  to  arc  spoking.  When  the  averaged  temperature  and  uncertainties 
from  Pulse  19  are  used,  the  mass  loss  drops  to  0.0176±0.046  pg.  This  suggests  that  in 
areas  where  the  arc  is  not  directly  over  the  propellant  face,  the  expected  amount  of  Teflon 
liberated  from  the  surface  is  extremely  low  for  a  given  pulse.  Therefore  the  total  mass 
likely  lost  after  the  pulse  is  given  an  upper  bound  of  1.2±0.6  pg. 
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5.  CONCLUSIONS 


This  paper  described  the  construction,  calibration,  and  application  of  HgCdTe 
photovoltaic  detector  diagnostic  to  determine  the  surface  temperature  of  Teflon 
propellant  during  and  after  the  pulse  of  a  microPPT.  This  diagnostic  is  shown  capable  of 
measuring  Teflon  surface  temperature  in  real-time  after  the  completion  of  the  current 
pulse.  A  theoretical  description  was  developed  to  predict  the  detector  voltage  output. 

This  theory  was  confirmed  by  predicting  experimental  calibration  data  and  then  used  to 
estimate  a  wavelength-dependent  emissivity  of  Teflon  in  the  LWIR.  A  calibration 
procedure  was  created  for  materials  with  wavelength  dependent  emissivity  and  high 
transmissive  character.  This  allows  application  of  emissive  temperature  measurements  to 
materials  that  until  now  have  been  difficult  to  characterize.  Using  data  from  infrared 
spectroscopy  and  infrared  thermography,  a  method  of  estimating  the  wavelength- 
dependent  emissivity  has  been  developed.  To  our  knowledge  these  are  the  first 
measurements  of  Teflon  surface  temperature  with  infrared  thennography  as  well  as  the 
first  application  of  high-speed  thennography  to  discharge  physics  for  analyzing  plasma- 
wall  interactions.  Although  data  during  the  pulse  of  the  micro-PPT  is  not  resolved  due  to 
plasma  interference,  the  cooling  curve  of  the  propellant  is  characterized  immediately  after 
the  pulse.  Measurement  of  the  surface  temperature  of  the  Teflon  propellant  during  the 
pulse  may  be  possible  if  the  line  transitions  in  the  plasma  are  sufficiently  different  from 
the  major  emissive  structures  in  the  solid  Teflon.  In  this  case,  a  line  filter  at  8.4  pm  over 
the  detector  may  remove  plasma  emission  sensed  through  the  broadband  nature  of  the 
detectors  used.  The  possible  limitation  expected  will  then  be  signal-to-noise  ratio. 
Predictions  of  surface  temperature  from  the  Keidar-Boyd  model  are  in  good  agreement 
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with  experiment.  Arc  constriction  modeling  shows  a  feasible  explanation  for  the  range  of 
surface  temperatures  measured  from  pulse-to-pulse  and  is  supported  by  DICAM  pictures 
of  anode  spotting  during  operation.  For  mass  loss  calculations,  an  assumption  of  full 
surface  coverage  is  used  to  provide  an  upper  bound  of  1.2±0.6  pg  which  corresponds  to 
23±1 1%  of  the  total  mass  lost. 
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Figure  Captions 


Figure  1 :  Linear  HgCdTe  detectors  shown  on  a  micro-PPT  face.  The  small  squares  are 
the  detector  images  and  the  gray  circles  are  the  micro-PPT  electrodes.  This  drawing  is  to 
scale. 

Figure  2:  For  calibration,  Teflon  is  sandwiched  between  a  heated  copper  plate  and  a  back 
plate.  There  is  no  heat  source  behind  the  optical  access  holes. 

Figure  3:  l/8th  inch  thick  Teflon  is  sandwiched  between  copper  plates  in  the  calibration 
jig.  Thermocouples  are  embedded  in  the  sides  to  approximately  the  optical  viewing 
location. 

Figure  4:  Schematic  of  the  (a.)  viewable  cone  for  the  Plume  measurement  technique  and 
(b.)  target  location  for  the  Surface  measurement  technique. 

Figure  5:  Estimated  a.)  non  plasma-exposed  Teflon  emissivity  and  b.)  plasma-exposed 
Teflon  emissivity  for  2-12  pm  wavelengths. 

Figure  6:  Calibrations  were  performed  with  virgin  Teflon  and  Teflon  exposed  to  micro- 
PPT  plasma  conditions.  The  plasma-exposed  data  is  used  for  calibration. 

Figure  7:  Comparison  of  surface  and  plume  measurements  in  a  micro-PPT.  4.35  J 

Figure  8:  A  test  with  14  pulses  shows  the  spread  of  cooling  temperatures  measured. 

Valid  data  begins  after  the  current  pulse  ends  (26  ps). 
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Figure  9:  Intensified  images  showing  arc  behavior  at  a.)  2.3  J,  b.)  5.55  J,  and  c.)  6.73  J. 
Capacitance  is  0.417  pF. 

Figure  10:  Measurement  of  a.)  the  highest  temperature  case  from  the  data  set  shown  in 
Figure  8  with  overlying  prediction  from  the  Keidar-Boyd  model  of  predicted  surface 
temperature  under  the  arc  constriction,  b.)  shows  the  lowest  cooling  temperature  from 
the  data  set  shown  in  Figure  8  with  overlying  prediction  or  surface  temperature  far  from 
the  arc  constrition  location. 

Figure  1 1 :  Surface  temperature  and  thruster  current  for  Pulse  20.  From  26-40  ps  the 
Teflon  vapor  pressure  is  calculated.  Data  before  20  ps  is  not  valid. 

Figure  12:  The  total  number  of  molecules  and  mass  evaporated  based  on  the  surface 
temperature  measurement  of  Pulse  20. 
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Figure  2:  Antonsen  et  al.  “Fast  Surface  Temperature. . 


Embedded  Thermocouples 


Figure  3:  Antonsen  et  al.  “Fast  Surface  Temperature. . 


a.  Plume  Measurement 


b.  Surface  Measurement 


Figure  4:  Antonsen  et  al.  “Fast  Surface  Temperature. . 
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Figure  5:  Antonsen  et  al.  “Fast  Surface  Temperature. . 
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Figure  6:  Antonsen  et  al.  “Fast  Surface  Temperature. . 
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Figure  7:  Antonsen  et  al.  “Fast  Surface  Temperature. . 
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Figure  8:  Antonsen  et  al.  “Fast  Surface  Temperature.. 
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Figure  9:  Antonsen  et  al.  “Fast  Surface  Temperature. . 
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Figure  10:  Antonsen  et  al.  “Fast  Surface  Temperature.. 
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Figure  1 1 :  Antonsen  et  al.  “Fast  Surface  Temperature. . 
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Figure  12:  Figure  10:  Antonsen  et  al.  “Fast  Surface  Temperature...” 
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